Identity method - a new tool for studying chemical fluctuations 
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Event-by-event fluctuations of the chemical composition of the hadronic system 
produced in nuclear collisions are believed to be sensitive to properties of the tran- 
sition between confined and deconfined strongly interacting matter. 
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fSj ■ In this paper a new technique for the study of chemical fluctuation, the identity 

method, is introduced and its features are discussed. The method is tested using 
data on central Pb-Pb collisions at 40A GeV registered by the NA49 experiment at 
<N i the CERN SPS. 
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I. INTRODUCTION 

The most interesting features of the phase diagram of strongly interacting matter are the 
Critical Point (CP) and the I s * order phase transition line. Event-by-event fluctuations of 
the chemical composition of the hadronic system produced in nuclear collisions are believed 
to be sensitive to both of them. The first relevant measurements were performed by the 
NA49 experiment at the CERN SPS. A systematic scan in beam energy and system size 
was recently started by the NA61 collaboration. Furthermore, additional insight is expected 
from the RHIC beam energy scan program. 

There are several measures used to quantify chemical fluctuations, among them: <J<i yn 



3| used by NA49 and v& yn J4] used by STAR. Both are related as o\ yn ~ fdyn and share 
the same disadvantage, namely, they depend on volume (number of wounded nucleons) 
and volume fluctuations in thermodynamical (wounded nucleon) models. Another measure 



6[ was used by NA49 to characterize transverse momentum (?], Q], electric charge 9| 



and azimuthal angle fluctuations [10|. It is free of the mentioned disadvantage of Gdyn, Vd yn - 
However, all these measures of chemical fluctuations are affected by non-perfect particle 
identification. This is illustrated in Fig. 1 (left panel), where the spectrum of the specific 
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energy loss (dE/dx) of particles measured by the NA49 Time Projection Chambers is 
shown for a selected phase-space bin. The dE/dx signal depends on the particle mass and 
together with the particle charge measurement is used to identify particles. It is seen that 
the dE/dx distributions of different particle species partly overlap and thus unique particle 
identification is not possible. 

The identity method adapts the <& x measure to take into account a non-unique particle 
identification, while keeping the advantages of <& x . 

II. IDENTITY METHOD 

Let us assume that particles are identified according to their measured mass. The mea- 
sured mass spectra of all particles and of particles of type h in the analyzed event sample 
are denoted as p and pt, respectively. The spectra are normalized to the corresponding 
mean multiplicities per event, namely: 

Jdmp(m) = (N) , Jdmp h (m) = (N h ) . (1) 

Furthermore, we define a single particle variable called the particle identity as: 

I s. def Phijn) 

Wh{m) = — — . 2 
p(m) 

The fluctuation measure *& w h is then introduced in a way similar (the roots over the two 
components are absent) to the $ measure. First, a single particle variable z is defined as: 

z = w h - wK , (3) 

where the bar denotes the inclusive mean and thus w = (N^) / (N) . Second, an event variable 
Z, which is the multiparticle analog of z, is calculated as: 

N 



Z = J2( W h( m i) ~Wh) , ( 4 ) 



def 

i=l 

where N is the multiplicity and i is the particle index in an event. 
Finally, the fluctuation measure ^f w h is defined as: 



, T , def (Z 2 ) -2 

y wh = -T^r - z 2 . (5) 



For further analysis one denotes two possible values of ^ w h- 
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• ^res which is the value of ^ w h for the experimental mass resolution case, 

• ^ corr which is the value of ^! w h for the perfect mass resolution case. 

In order to correct for the non-unique particle identification we calculate the variance per 
particle due to random identification for the experimental mass resolution case: 

1 f°° 

Var res = — — / dmp(m) ■ w h (m)(l - w h {m)) . (6) 
\JS ) Jo 

It is easy to show that for unique particle identification (the perfect mass resolution case, 
Whijn) = 8{m — nth)) the result is Var res = f Var^ = 0, whereas for no mass resolution ( 
Wh(m) = const. ) one obtains Var res = f Var B = ^^-(1 — (Nf l )/(N)). For the experimental 
data analysis the integral in Eq. [6] is replaced by a sum over all particles. 
The following key relation can be proven 1 : 

^corr = ^res(l - V ar res /V ar B y 2 . (7) 

It shows that the measured fluctuations ^ res can be corrected for the effect of non-unique 
particle identification in a model independent way. This is because the correction factor, 
(1 — Var Tes /VarB)~ 2 1 depends only on the experimental resolution and mean particle mul- 
tiplicities. 

Equation [7] was checked by numerous Monte Carlo simulations with different types of 
correlations, mass resolution functions and mean particle multiplicities. The results of these 
simulations are shown in Fig. [1] (right panel). 



III. THE IDENTITY METHOD TEST USING NA49 DATA 

In the analysis of experimental data we use the particle energy loss dE/dx as the measure 
of the mass m. For optimal identification the dE/dx spectra from the NA49 TPCs are 



determined in bins of tota 
electric charges separately 



and transverse momentum, azimuthal angle as well as for both 

1 



-|3J. In each bin four Gauss functions (for electrons, pions, kaons 
and protons) are fitted. An example of such a fit is displayed in Fig. 1 (left panel). The 
fitted functions are then used as the ph and p functions of the identity method. The further 
analysis steps are as follows: 
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using mean particle multiplicities the variance Var B = (N h )/(N)(1 — (N h ) / (N)) is 
obtained, 

for each particle its identity is calculated 

. | ,x p fe (< rfff/cb >j \q,Ptot,PT,<P) , Q s 

w hi (< dE/dx>i\q,ptot,PT,(p) = / at? I a ^ I TT 5 ( 8 ) 

p(< dE/cb >j \q,Ptot,PT,<P) 

using the experimental dE/dx resolution functions, p/j and p (M is the total number 
of particles used in the analysis) the variance 

M 

M 



Var res = -jr^WhiiK dE/dx > { \q,p to t,PT,(f>)(l -%(< dE/dx >» |g,Ptot,Pr») 



l 



(9) 

is computed, 

• using the identity values, Whi, ^res is calculated and 

• the corrected value of ^ corr is obtained using Eq. [3 

As a first test of the identity method proton fluctuations were studied in Pb+Pb collisions 
at 40A GeV energy. Positively and negatively charged particles with total momentum up 
to 40 GeV/c and transverse momentum up to 2 GeV/c were used for the analysis. The 
total number of analyzed events was 4000. The mean multiplicities are < N >= 165.40, 
and < N p >= 42.16. The obtained value of \& res ■ 1000 = —17 A ± 3.5. The corresponding 
correction factor for non-unique particle identification was calculated to be ~ 1.2. 
The value corrected for the finite resolution is \l/ corr ■ 1000 = —22.3 ± 4.4. The analysis of 
Pb+Pb collisions at all NA49 energies is in progress. 
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FIG. 1: Left panel: Distribution of specific energy loss measured in the NA49 TPCs for positively 
charged particles is a bin p tot e(4.4 — 5.3) GeV/c, pTe(0.0 — 0.2) GeV/c and (/>e(0.757r — ir). The fitted 
Gauss functions are shown by solid curves. Right panel: The ratio ^res/^corr 

versus Var res /VarB 

calculated within several Monte Carlo simulations with different parameters of experimental res- 
olution, particle multiplicities and fluctuations. The results agree with the analytical dependence 
given by Eq. [7J 
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FIGURE CAPTIONS 

1. Left panel: Distribution of specific energy loss measured in the NA49 TPCs for pos- 
itively charged particles is a bin p tot e(4.4 — 5.3) GeV/c, pre(0.0 — 0.2) GeV/c and 
0e(O.757r — it). The fitted Gauss functions are shown by solid curves. Right panel: 
The ratio ^ res /^ corr versus Var res /VarB calculated within several Monte Carlo simu- 
lations with different parameters of experimental resolution, particle multiplicities and 
fluctuations. The results agree with the analytical dependence given by Eq. 



